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Abstract 

A continuum model of crystalline solid equilibrium is presented in 
which the underlying periodic lattice structure is taken explicitly into 
account. This model also allows for both point and line defects in the 
bulk of the lattice and at interfaces, and the kinematics of such defects 
is discussed in some detail. A Gibbsian variational argument is used 
to derive the necessary bulk and interfacial conditions for multi-phase 
equilibrium (crystal-crystal and crystal-melt) where the allowed lattice 
variations involve the creation and transport of defects in the bulk and 
at the phase interface. An interfacial energy, assumed to depend on the 
interfacial dislocation density and the orientation of the interface with 
respect to the lattices of both phases, is also included in the analysis. 
Previous equilibrium results based on nonlinear elastic models for inco- 
herent and coherent interfaces are recovered as special cases for when the 
lattice distortion is constrained to coincide with the macroscopic defor- 
mation gradient, thereby excluding bulk dislocations. The formulation is 
purely spatial and needs no recourse to a fixed reference configuration or 
an elastic-plastic decomposition of the strain. Such a decomposition can 
be introduced however through an incremental elastic deformation super- 
posed onto an already dislocated state, but leads to additional equilibrium 
conditions. The presentation emphasizes the role of configurational forces 
as they provide a natural framework for the description and interpretation 
of singularities and phase transitions. 

1 Introduction 

Consider two phases of a given material separated by a sharp interface and in 
equilibrium in a heat bath. Let one phase be a crystalline solid, the other phase 
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may be either the corresponding fluid melt or another crystalline solid (such 
as in a twinned crystal). This paper is concerned with the determination of 
the equilibrium conditions for the bulk regions occupied by the two phases and 
for the interface separating them — a problem of significant interest in diverse 
fields. Numerous authors have already studied various aspects of it and have 
derived equations governing the mechanical, thermal, and chemical equilibrium 
(Eshelby, 1970; Robin, 1974; Larche & Cahn, 1978,1985; Cahn, 1980, Cahn 
& Larche, 1982; Grinfel'd, 1981; Mullins, 1984; Alexander & Johnson, 1985; 
Johnson & Alexander, 1986; Kaganova & Roitburd, 1988; Leo & Sekerka, 1989; 
Gurtin, 1993; Cermclli & Gurtin, 1994; Leo & Hu, 1995). A fundamental 
characteristic of a crystalline solid is its lattice structure, in which the atoms are 
periodically spaced. Additionally the lattice contains defects, e.g., point defects 
(interstitials and vacancies) and line defects (dislocations), both types of which 
can impact on the equilibrium behavior. Larche & Cahn (1978) included point 
defects in their original analysis by treating the interstitials and vacancies as 
additional mobile species subject to a lattice or network constraint. This work 
was later extended to allow for surface energy associated with the interface 
(Cahn, 1980, Cahn & Larche, 1982; Mullins, 1984; Alexander & Johnson, 1985; 
Johnson & Alexander, 1986; Leo & Sekerka, 1989; Gurtin, 1993; Cermelli & 
Gurtin, 1994; Leo & Hu, 1995). In this approach however the lattice does not 
appear explicitly, so it is not clear how other properties of the lattice, such as the 
distortion of the lattice or the occurrence of dislocations, affect the governing 
equations. 

In this paper we present a continuum model of crystalline solid equilibrium 
that explicitly incorporates the underlying lattice structure and includes the 
possibility of both point and line defects. Following Larche & Cahn (1978), we 
then give for this model a Gibbsian variational analysis, in which we obtain 
the necessary conditions for two-phase crystal-crystal equilibrium. We consider 
general non-equilibrium variations that allow for creation and transport of de- 
fects in the bulk and at the interface. We treat both coherent and incoherent 
interfaces, two extreme cases that characterize many types of interfaces between 
crystalline solids. In particular, incoherency is important since it can arise from 
the presence of defects at the phase interface. Additionally, we consider the 
case of crystal-melt equilibrium. We also include interfacial energy, assuming 
that it depends on the interfacial dislocation density and the orientation of the 
interface with respect to the lattices of both phases. 

We model a crystalline solid by a continuum microstructural theory where 
the microstructural variables represent in an average sense the underlying lattice 
structure: 

- 3 linearly independent vector fields denote the average lattice vectors at a 
spatial point; 

- interstitials and vacancies are treated in the usual way by modelling them 
as additional mobile species subject to a lattice constraint; 

- dislocations are related to the anholonomicity (i.e., non-integrability) of 
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the micro-structural lattice vectors. 



In this model we introduce neither a fixed reference configuration nor a macro- 
scopic deformation gradient; rather we describe the state only by the microstruc- 
tural variables through a spatial description, where quantities are defined on a 
per unit cell basis at a fixed point in space. A spatial description is uncommon 
in traditional models of solids, but is more convenient when interfaces are in- 
coherent as they remain together in the current configuration. It also seems to 
be the most appropriate when relating to standard x-ray observations of lattice 
points, which do not indicate how the lattice points actually evolve. We also 
assume that the energy per unit cell of the current configuration depends upon 
the local current state (i.e., lattice vectors, dislocation density, and occupancy 
density). 

Commonly in microstructural models, lattice vectors are assumed to deform 
as material line elements — the so-called Born rule (Ericksen, 1984) — so that the 
lattice vectors in the deformed configuration are related to the lattice vectors 
in the reference configuration by the macroscopic deformation gradient. This 
assumption however is not always valid: in particular it fails when there are 
dislocations; it can also fail in other cases (Zanzotto, 1992). We therefore for- 
mulate our model without employing the Born rule, but then we show how 
more classical models for solids are recovered when the Born rule is imposed as 
a constraint. Thus our equilibrium results are new in that 

(1) bulk and surface dislocations are explicitly taken into account; and 

(2) no a priori relation between reference lattice vectors and actual lattice 
vectors is assumed. 

When the Born rule is imposed as a constraint, the lattice vectors then deform 
as material line elements, the dislocation density vanishes in the bulk, the lat- 
tice distortion becomes equivalent to the deformation gradient, and we recover 
equilibrium results equivalent to those of Larche & Cahn (1978), Leo & Sekerka 
(1989), and Ccrmclli & Gurtin (1994). 

Our approach, however, does not require any recourse to the concept of 
macroscopic deformation or fixed reference configuration. On the other hand, 
for a number of applications or in certain special cases, it may be convenient to 
introduce the often used elastic-plastic decomposition of the strain (Bilby et al, 
1957; Kroner, 1960; Lee, 1969). We do so by identifying the dislocated crystal 
with the intermediate configuration of such a decomposition and then superpos- 
ing a classical elastic (i.e., defect-preserving) deformation onto this dislocated 
state. In this approach our "elastic" part of the elastic-plastic decomposition is 
truly elastic in the traditional sense. This procedure allows one to construct a 
simple model of dislocated bodies while keeping track of all the different elastic 
and plastic components of the theory. Due to the additional structure in the 
theory arising from the incremental elastic deformation, additional equilibrium 
conditions arise, so that the resulting equilibrium theory is not equivalent to that 



3 



without the elastic-plastic decomposition. In particular, when such an interme- 
diate configuration is introduced, both the dislocated intermediate configuration 
and the current configuration must be in mechanical equilibrium. 

There is some precedence for a microstructural approach to modelling crys- 
talline solids. Motivated by the Cosserat director theories for oriented materials 
(Ericksen & Truesdell, 1958; Truesdell & Toupin, 1960), both Fox (1966,1970) 
and Toupin (1968) proposed using a triad of continuous vector fields to rep- 
resent lattice vectors associated with a material point and demonstrated how 
dislocations could be related to the non-integrability of the directors. The de- 
formation associated with the lattice directors was assumed distinct from the 
average, macroscopic deformation. Additionally Fox discussed how a notion of 
lattice slip could be introduced. Mullins & Sekerka (1985) also used a triad of 
vector fields to model the lattice vectors, which, however, where taken to de- 
form as material line elements (the Born rule). Consequently any constitutive 
dependence on them could be replaced by a dependence on the macroscopic de- 
formation gradient, which excluded the possibility of bulk dislocations. In this 
case their model described anisotropic elastic solids (Ericksen & Rivlin, 1954). 
Davini (1986) and Davini & Parry (1989,1991) discussed in detail the various 
defect measures and their relation to the lattice directors. They also introduced 
the neutral deformations — inelastic deformations that do not change the de- 
fect content — and demonstrated a corresponding decomposition of the strain. 
Naghdi & Srinivasa (1993a, b, 1994a, b) formulated a director theory as a model 
of slip in plastic deformations. Similarly, Bcsscling and van der Giessen (1994) 
introduced a triad of vector fields for the lattice vectors in their discussion of 
inelastic deformations, but did not introduce the notion of defects in this model. 
More recently, Dluzewski (1996) discussed the driving forces acting on defects in 
a director model. And Ericksen (1997) discussed an equilibrium theory of crys- 
tals in terms of lattice directors without employing the Born rule and without 
introducing a macroscopic deformation gradient, but did not consider defects. 
Our approach follows closely that of Davini, Parry, and Ericksen in that it is 
based entirely on a description in terms of the underlying lattice structure. 

Finally, our presentation emphasizes the role of the configurational forces^ as 
they provide a natural framework for studying singularities and phase transitions 
(Eshelby, 1951,1970; Maugin, 1993; Gurtin, 1995). Configurational forces do 
work over changes in the reference lattice and arise naturally in our theory 
as the derivatives of the energy with respect to the reciprocal lattice vectors 
(Cermelli & Sellers, 1998). In our calculations the configurational forces result 
as primary quantities and provide an immediate interpretation of the resulting 
equilibrium conditions. 

2 Crystalline Solid Model 



-"^Also called Eshelby forces or energy-momentum. 
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2.1 Motivation 



By the expression crystalline solid, we mean a model for a solid that explicitly 
incorporates the underlying lattice structure typical of a crystal. Our model of 
a defective crystalline solid is based on the notion that the regular lattice struc- 
ture is, on the average, locally recognizable and measurable (Davini, 1986). For 
instance, x-ray observations can provide local values for the various lattice pa- 
rameters. Further, such experimental measurements indicate that the dilatation 
of the lattice parameters does not always correspond to the average macroscopic 
strain as measured in, for example, tensile tests. These results intimate that 
the macroscopic deformation gradient can be in some cases an inappropriate 
measure of the true deformation of the crystal; finer kinematic details of the 
lattice may be needed. In the following we introduce microstructural variables 
to represent such kinematic details of the lattice and its distortion. In this point 
of view, the state of a crystalline solid with defects is completely determined 
by the crystalline structure, which we take as the set of lattice vectors and the 
occupancy density of lattice sites. These quantities are viewed as fields over the 
region of space occupied by the material in its actual or current configuration. 

2.2 Lattice microstructure 

Let Scl 3 denote the region of space occupied by the crystalline solid in its 
actual or current configuration. We will consistently use this spatial region as 
a reference. 

To each spatial point x 6 £>, we associate a triad of vector fields e^(x) (where 
i ranges from 1 to 3) that represent the lattice basis vectors at x. They describe 
in an average sense the local microscopic arrangement of the lattice points and 
vary continuously from one spatial point to another. It is also useful to introduce 
the reciprocal lattice vectors e l (x) defined by the relation 

e* • Bj = 8$ , e» <g> e' = e' ® = 1, (1) 

where <5* is the Kronecker delta symbol and repeated Latin indices are summed. 
Generally x-ray observations provide direct measurements of the reciprocal lat- 
tice vectors e\ We shall also denote a fixed, constant and uniform lattice basis 
by Ei, with reciprocal basis E\ It provides a convenient reference lattice from 
which to measure the change of the lattice vectors.^ 
The lattice distortion is by definition the tensor field 

Fi at =e 4 ®E l , (2) 

which takes the reference lattice vectors Ej to the actual lattice vectors ej at 
each spatial point. Similarly, the inverse lattice distortion 

Ff^E^e* (3) 

2 Although we do not introduce a fixed reference configuration for the current configuration 
B, we do use the notion of reference lattice. 
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takes the actual lattice vectors to the reference lattice vectors at each spatial 
point. We do not assume that Fi at or Fj~J are gradients. In particular, Fi at 
need not correspond to the macroscopic deformation gradient F. 

Since the reference lattice vectors are constant, the actual or current state 
of the lattice is equivalently characterized by either the fields e,*, here taken as 
spatial fields on £>, or by the fields e 1 , due to the relation ([[]) between lattice 
vectors and reciprocal lattice vectors. 

Furthermore, we interpret ei • (e2 x 63) as the volume of a unit cell, so that 

n := e 1 ■ (e 2 x e 3 ) = — ^ (4) 

ei • (e 2 x e 3 ) 

is the number of cells per unit volume. It is also a continuously varying spatial 
field defined on B. 



2.3 Point defects 

In a perfect crystal the atoms occupy lattice points. In a defective crystal 
some of the lattice points may not be occupied, these being called vacancies. 
Additionally, some atoms may be at points not located at the lattice points, 
these being called inter stitials. Here, to maintain the presentation simple, we 
assume that there exists a single mobile species of atoms, either vacancies or 
interstitials. The extension to an arbitrary number of species is straightforward. 
We thus introduce the following scalar fields: 

g(x) atoms per unit cell 

£Vac(x) vacancies per unit cell 

If Qv&c > 0, we will interpret g vac as a density of vacancies; if g vac < 0, we 
will interpret |g V ac| as a density of interstitial atoms. Since we allow for only a 
single mobile species, only the sign of g vac distinguishes between vacancies and 
interstitials. These two fields are not independent but are related by the lattice 
constraint: 



g + gvac =£ = const., (5) 

where £ is the number of lattice points per unit cell. The lattice points are 
regarded as a geometric constant of the lattice. Consequently only one of the 
quantities is independent, and we arbitrarily choose g as the independent den- 
sity. 

If independent mobile point defects are allowed, they must also be specified 
in addition to the lattice vectors. Thus for our model, the current state of 
a crystalline solid with a single mobile species is given by (e^p), specified as 
spatial fields on B. Alternatively and equivalently, the lattice vectors ej can be 
replaced by the reciprocal lattice vectors e l or the molar density g replaced by 
the vacancy density g V ac- 
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2.4 Line defects 



Classically, the Burgers vector b = b % Ej is taken as a measure of the dislocations 
in the sense thatQ 

b*= I e* dl= /(curl e>n c cfo, (6) 

represents the atomic lattice displacements in the reference lattice for a circuit 
dC of a fixed surface CcB with unit normal n c (Bilby & Smith, 1956; Davini, 
1986). To see this, notice that the quantity b 1 defined by (0) is the net number 
of i-planes (i.e., planes orthogonal to e*) which are crossed by the curve dC, 
so that any term in the combination b = fr J Ej represents: (the net number 
of added i-planes encountered upon traveling along dC) times (the reference 
i-lattice vector) = (the net Burger vector of the circuit). 

By definition the reference lattice vectors Ej are constant and dislocation 
free. Define the 3 quantities 

g*(x) := curl e*(x). (7) 

Commonly Ej ® g l is used as a tensorial measure of the density of dislocations 
and corresponds to the dislocation density introduced by Nye (1953). Since 
the reference lattice vectors Ej are constant, we use instead the 3 vectorial 
dislocation densities g l , which are taken as a spatial fields on B. 
We will frequently use the equivalent fields 

«^:=I g W, (8) 
n 

which can be interpreted as the dislocation density per unit cell (Davini, 1986; 
Davini & Parry, 1991). 

We will model phase boundaries as sharp surfaces. To allow for this possi- 
bility, suppose now that the circuit dC is intersected by a sharp surface S across 
which there is a jump in lattice vectors. Then (J^j must be modified to 



b l = / e 1 -d\ = / (curl e 4 )-n c da- / [e l ] • dl, (9) 
Jdc Jc Jdcns 

where J-J denotes the jump of a quantity across the surface S in which the jump 
is taken as the limit from the portion of B into which the normal n s to the 
surface points. Since the above relation is independent of the surface C (keeping 
fixed the boundary curve dC), the tangential jump in the lattice vectors provides 
a measure of the defectivity of the interface. Thus we introduce 

[Pel (10) 
as 3 vectorial measures of the surface dislocation density, where 

P:=l-n s ®n s (11) 



3 In accord with our use of the spatial description, the derivatives are with respect to the 
current position x. 
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is the projection operator onto the tangent plane of S. The quantity [E; <g) 
(Pe 2 )]] corresponds to the tensorial surface density introduced by Bilby (1955). 
It is equivalent to (|Io| ) since the reference lattice vectors are continuous and 



linearly independent. To see that (10) indeed measures the density of interfacial 



dislocations, let u be a tangent vector to the interface S. Then 

[E i ®(Pe i )lu=([e i ]-u)E i . (12) 

And since (e*) a - u and (e 1 )/ 3 ■ u are the net number of i-planes from side a 
or (3 which intersect the vector u (recall that |e l | is the reciprocal of the dis- 
tance between two consecutive i-planes), [e 4 ]] • u is the net number of i-planes 
intersecting u. Thus the sum ( |l2| ) consists of the net number of added i-planes 
encountered upon traveling along u times the reference i-lattice vectors, which 
equals the net Burgers vector of the interfacial dislocations along u. 

In classical elastostatics, a deformation is holonomic (i.e., compatible or 
integrable) if the strain tensor is the gradient of a global, smooth map. For our 
model of a crystalline solid, this notion is expressed by a compatibility condition 
stating that there exists a smooth map x such that 

Ei g> e* = grad X - (13) 

In other words, the state of a crystalline solid is holonomic if the inverse lattice 
distortion Fi at is the gradient of a global, smooth map. In this case, the notion 
of deformation arises naturally since x corresponds to the macroscopic inverse 
deformation to the reference lattice. We will call a lattice a holonomic lattice 
(with respect to E^) if ([l3|) is satisfied for some X- The reference lattice, by def- 
inition, is holonomic; the distorted lattice however is in general not holonomic. 
And since the reference lattice vectors Ej are constant, (|l3| ) shows that the bulk 
dislocation densities g l vanish identically for a holonomic lattice (though the 
surface dislocation densities ( |l0| ) need not vanish at a singular surface S of a 
holonomic lattice). In this model of a crystalline solid, dislocations are viewed 
as the obstruction to patching together local measurements of the reciprocal 
lattice to a single, global lattice. 



2.5 Born rule 

The Born rule states that lattice vectors deform as material line elements under 
a given deformation (Ericksen, 1984). We are primarily concerned with the case 
when this does not hold, as the Born rule excludes dislocations in the bulk. We 
will consider however the consequences of imposing it as a constraint, where it 
will be shown to imply common models for elastic equilibrium as special cases. 

Consider now a holonomic lattice in the actual configuration, so that the 
dislocation density vanishes inside B. There is a natural inverse deformation % 
for the lattice, through which we can associate %(6) with the region occupied 
in the undeformed configuration. In this case, the Born rule yields (cf. jl^)) 

ei = (grad x _1 )Ei or e* = (grad X ) T E\ (14) 
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If we assign on x{&) the perfect lattice determined by E 1 , then (|l4|a) states that 
the lattice vectors deform as material vectors through Alternatively, (|l4|b) 
states that the reciprocal lattice vectors deform as material vectors through X- 
Using the Born rule at a surface of discontinuity S, the 3 surface dislocation 
densities (|Io| ) become 

[P(grad X ) T flE l , (15) 

with P the projection operator onto the tangent plane of S. They correspond 
to the single tensorial quantity [P(grad x) T ] used by Cermelli & Gurtin (1994). 



3 Lattice Variations 

Let (e i; g) denote as before a state of a body occupying the region B. We 
denote an arbitrary smooth variation of this state by (Sei,5g), which does not 
necessarily correspond to any equilibrium state.^ The external boundary dB is 
assumed rigid and fixed, so the variations vanish on dB (but not necessarily on 
an internal boundary such as the phase interface). Consider now the additive 
composition of the two configurations 

e lX := e { + XSe^ g x :=g + XSg (16) 

where A is a small parameter. This composition induces a corresponding varia- 
tion in any functional $ depending on the composed configuration: 

(17) 

In particular, it induces variations in dislocation density g 1 and the reciprocal 
lattice vectors e 4 , related by the identity 

(5g 4 = curl 5e\ (18) 

The induced variation of the Burgers vector is Sh = Sb l Ej where 

Sb*= [ Sj-dl. (19) 

J dC 

Generally the Se l induce a dislocation flux (Kosevich, 1962). 
A straightforward calculation yields the relations 

Sn = nei • Se l , 

n 5g i3 = - (g< ■ e j )e ; • Se l + e j ■ 5g l + g l • Se j , (20) 
which will be needed later. 

4 Since we are using a spatial description, the variations correspond to the Eulerian varia- 
tions of Leo & Sekerka (1989) 



8® = ^*[ e *A, Q\] 
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Consider now an internal boundary given by the singular surface S. We allow 
the e, and the dislocation density to be discontinuous across this surface. A 
variation of the position of the interface in space is described by 

5r = n s 5r, (21) 

with n s a choice of unit normal to S. Variations of the position of the surface 
S and of lattice vectors at the surface will involve in general the creation and 
transport of defects. 

4 Equilibrium Equations without Interfacial En- 
ergy 

In section 2 we developed the kinematics to describe an equilibrium state of 
a crystalline solid with point and line defects. In section 3 we extended the 
kinematic structure in order to discuss variations of states that induce creation 
and transport of such defects. In this section we employ this kinematic struc- 
ture to derive the equations governing multi-phase equilibria. The approach is 
variational and involves rendering stationary an appropriate energy potential. 
We treat crystal-crystal and crystal-melt equilibria separately. Also we consider 
first the case without interfacial energy then in the following section extend the 
equilibrium results to include interfacial energy. 

4.1 Crystal-crystal equilibrium 

Let two crystalline solid phases (denoted by a and (3) share a common interface 

5 and be immersed in a heat bath at temperature 9b- Further, let the external 
boundaries be impermeable and fixed. Such a situation can be in equilibrium 
only if an appropriate thermodynamic potential is stationary. We now construct 
such a potential and determine the necessary conditions for it to be stationary. 
We continue to use a spatial description and take all extensive quantities de- 
scribing the crystalline solid as per unit cell in the actual state. 

Since there is no external mass supply the total mass of the system is con- 
stant and is given by 

M = M a + M = { ngdv. (22) 

JB a VB^ 

We assume that the total internal energy of the system is given by 

E = E a +E< 3 = f nedv, (23) 

where e is the internal energy density per unit cell in the actual state, and any 
surface energy has been assumed negligible. Furthermore, we assume that the 
entropy is additive: 

S = S a + S = f nsdv, (24) 



10 



where s is the entropy density per unit cell. 

A constitutive relation must specify one these quantities as a function of the 
current state of the material. As a local measure of the lattice structure, we 
choose the quantities (ei,g y , g), since ej specifies the periodic structure of the 
lattice, <? u the dislocation density per unit cell, and g the atoms per unit cell. 
This choice is not unique; it is also equivalent to (e 1 , g l , g vac ). As a measure 
of the local thermodynamic state, we choose the entropy s per unit cell. Our 
constitutive assumptions are then that the internal energy per unit cell of each 
phase is given by: 

e = e a (ei,g iJ , g, s) for phase a, 

e = e' 3 (e i , g lj , g, s) for phase (3. (25) 

Notice that we assume that the energy depends on the current state as de- 
termined by the local lattice structure and entropy. In particular, there is no 
dependence on a macroscopic deformation gradient since we have not introduced 
a fixed reference configuration for B a U fi". 

Following Gibbs (1878), the grand canonical potential 

Q:= E-6 B S - fiM (26) 

(at constant total mass and total entropy) is stationary in equilibrium, where 
fi is a constant Lagrange multiplier ensuring conservation of mass. A necessary 
condition for its stationarity is that the first variation vanish: 

5n = SE-e B 5S-fiSM = 0. (27) 

Thus using (|22|)-(p7|) along with ( pp| ) and denoting by 

the jump of a field h at the interface from the a and /3-sidesJ^] we have 



5Q = 5 / nedv — 9sS nsdv — fiS ngdv 

'e^UB" JB a UBP Jl3 a UBP 

de 

[n(e — Obs — fig) • 6e l + n(— fi) 5g 

8°u8' 3 "g 

/<9e „ s r de de , 

+ ~ B b)5s + n— ■ 6e t + n-— 5g J \ dv 

- I fn(e - 9 B s - fig)} Sr da = (28) 
Js 

where the surface integral arises from the variation of the interface and reflects 
the induced change in phase, i.e., accretion. 

The main difficulty is to choose an appropriate set of independent variations 
in the bulk and at the interface that yields equations in a physically transparent 
form. To this end, let 

lu := e — Obs — fig. (29) 



5 The choice of n s pointing outward from phase a is consistent with our previous convention 
on the use of H 
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denote the grand canonical potential per unit cell. Further we introduce 



ki~-^je j , t*:=n|p t< := (t j - ei)e 3 -, (30) 

and we identify E'Cgik^ as the dislocation couple tensor and e*(S>t, as the Cauchy 
stress tensor (Cermelli & Sellers, 1998). The choice of basis indicates that the 
Cauchy stress acts on the actual lattice, whereas the dislocation couple acts on 
the reference lattice. We also introduce the configurational stress tensor E l ®Ci 
whereQ 

d(nu>) 



de 1 



ncuei - tj - (g J x ej) x kj. (31) 



is a generalized Eshelby relation containing an additional term due to the dis- 
locations. Again the choice of basis indicates that the configurational stress, 
in contrast to the Cauchy stress, acts on the reference lattice. In terms of the 
quantities (^0|) and ( |3l| ) the variation (pffi) can be expressed alternatively as 

C da da 

SSI — / \n(- u)Sg + n(- 6b) Ss + (c{ + curl ki) ■ Se 1 ] dv 

J B c UB p dg ds 

- [ ([nw] Sr + n s ■ ISe 1 x k,]) da = 0, (32) 
Js 



where we have used the divergence theorem and (|l8j) . The expression (|32|) of the 
first variation of the energy functional is now in a convenient form for deriving 
the appropriate necessary conditions. 

The choice of the possible independent variations should reflect the rele- 
vant physical situation. We consider two different cases: one corresponding 
to no constraints on the possible variations, the other corresponding to when 
holonomicity is imposed as a constraint, thus excluding variations that create 
dislocations in the bulk. 

4.1.1 Case I: unconstrained variations 

• Bulk Conditions: We take 

8g, ds, 5e l (33) 

as the admissible class of independent variations in the continuous bulk regions. 
These variations allow defect creation and transport in the bulk. With ([52j) 
they yield the following necessary conditions: 



uniform chemical potential 

de a de 13 



de -8 Q - ^ (34) 



6 The generalized Eshelby relation ( pl| ) has been derived by an invariance argument in 
Cermelli & Sellers (1998). 
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uniform temperature 

-d7 = ^7 = dB > (35) 

dislocation couple balance for each phase 

Ci + curlkj=0. (36) 



The uniformity of the chemical potential ( p4| ) and temperature ( |3q ) are standard 
and represent chemical and thermal equilibrium. The dislocation couple balance 
( p6| ) represents mechanical equilibrium. Additionally, it implies immediately the 
configurational force balance for each phase 

div a = 0, (37) 

which, with the Eshelby relation ( |3l| ) , is itself equivalent to the Cauchy force 
balance 

div (e l ®ti) = 0, (38) 

which is the more standard mechanical equilibrium condition. But provided 
that the lattice variations Se 1 are arbitrary, relations (|3^)-(^6|) are the indepen- 
dent necessary conditions for equilibrium in the bulk. 

• Interfacial Conditions: We want to allow for defect creation and transport at 
the interface so that we take the variations of the lattice vectors Se 1 on each side 
of the phase interface as independent, but require that Sr be continuous across 
the interface since the two phases remain in contact. Thus the independent 
variations at S are 

Sr, (6e*) a , (^f- (39) 

They correspond to an incoherent interface. With this choice of independent 
variations at the surface, ( |32| ) yields as necessary interfacial conditions: 

- continuity of the grand canonical potential 

H = 0; (40) 

- vanishing of the individual tangential couples 

P(ki) a = and P(k l ) /3 = 0. (41) 



The continuity of the grand canonical potential fl40| ) is standard and is due to 
the fact that the phases may not separate at the interface. The vanishing of 
the individual tangential couples reflects the fact that an incoherent interface 
must be in mechanical equilibrium with respect to both lattices. Moreover, this 
latter condition implies in turn the 
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- vanishing of the individual conjigurational tractions 

(c 4 ) Q • n s = and (c t f ■ n s = 0. (42) 

As in the corresponding bulk conditions on the configurational stress, (^2|) are 
not independent conditions since the identity 



/ (c 4 ) tt • n s da = - f (k t ) a ■ dl, 

JS' JdS' 



(43) 



(which holds for any subsurface S' C S) and (|4l]) imply ( |42| ) . Additionally, the 
identity 

n s • le l ® ti]n s = [nw] - n s • [e* g> c,]]n s - n s • Je 4 <g> ((g 3 x e 4 ) x kj)]n s , 

(44) 

which follows from the generalized Eshelby relation |3l]), shows that the conti- 
nuity of the grand canonical potential is equivalent with ( fill) an d (^) to the 

- continuity of the normal Cauchy traction 

n s -Ie 4 ®t 4 K =0. (45) 

Thus the independent equilibrium conditions at the incoherent interface are ( (4l| ) 
and either (|4(]) or, equivalently, (|H|). 



4.1.2 Case II: holonomic variations (Born rule) 

Since dislocations are in general created or transported when the lattice varia- 
tions Sei are arbitrary, the configurational force balance (j3?j) and corresponding 
interface conditions (f|2]) are consequences of the dislocation couple balance ( |36| ) 
and corresponding interface condition (|4l|). We now introduce a restricted class 
of lattice variations that yield (|37| ) and ( p2|) as independent necessary conditions. 

Consider the class of holonomically constrained variations for which the Born 
rule ( |l4| ) holds: variations of the lattice basis Se l are such that curl (i5e 4 ) = 0. 
Thus away from the surface S the variations can be expressed in terms of the 
inverse deformation x by the relation 

5e l = (grad S X ) T E\ (46) 

Since the external boundary dB is fixed, 8\ vanishes at dB, which excludes any 
rigid body displacements. Thus the actual and varied states may be realized as 
deformations of the same perfect-lattice reference configuration (%(£?), Ei). In 
the terminology of Davini and Parry (1989), the actual and varied states are 
elastically related. This is false for arbitrary variations not vanishing at dB, 
since in that case rearrangements^ such as slip may occur, so that the states e* 
and e* + <5e* are only neutrally related. 

7 Rearrangements are changes of state which do not modify the local lattice microstructure. 
A discussion of the equilibrium conditions for a crystalline body when global rearrangements 
are allowed, is due to Fonseca and Parry (1992). 
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The Born rule excludes dislocations in the bulk but not necessarily at an 
incoherent interface. Remember that at a coherent interface, x 1S continuous, 
whereas at an incoherent interface, x is discontinuous. Thus with the Born rule 
(14), the continuity of x implies the vanishing of the surface dislocation density 
(10) whereas discontinuities in x ma Y lead to nonvanishing surface dislocations. 

In order to represent these two possibilities more succinctly, we introduce 
the quantities 

3sX l : = ($X + (g ra d ns x) 8r) ■ E*, grad ns x := (grad x)n s , (47) 

which are, respectively, the components of the variation of the inverse defor- 
mation x following the variation of the interface S and the derivative of the 
deformation normal to the current interface. Note that 5sX l is continuous when 
X is continuous: 

M = fe1=0. (48) 

By a procedure analogous as before, we may write the first variation of the 
energy functional in the form 

Sn= [n(— - fi)Sg + n(— - 9 B ) 6s + (div a) Sx • E l l dv 

Jbcubp oq OS 

- {{a S sx l l ■ n s - n s ■ [e* ® t^]n 5 5r) da = 0, (49) 



where now the configurational stress component is defined as 

d(nu>) 

so that 



(50) 



(grad x) T (E i ®c l ) ^nwl-e^ti (51) 

corresponds to the usual Eshelby relation as expressed in the spatial description. 
As there are now no dislocations in the bulk, there is no dislocation contribution 
to (|5l]) and the configurational stress coincides with the classical Eshelby stress. 

• Bulk Conditions: Taking as independent variations in the bulk 



Sq, Ss, Sx, (52) 
we obtain now the following necessary conditions: 



uniform chemical potential (34); 
uniform temperature (|35|); 

configurational force balance (^) for each phase. 
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Because of the Eshelby relation (|5^) , the configurational force balance (^) in 
the bulk is equivalent to the Cauchy force balance (]3^). But for this holonomic 
case, the dislocation couples k; are constitutively indeterminate and are now 
defined by the couple balance Q36 ) . Thus the independent conditions with the 
class of holonomic variations are ( 34|) , (^) , and either (^7|) or ([38]) . 

• Incoherent interfacial conditions: We want to allow for surface dislocations 
at an incoherent interface, so that we take as the admissible class of independent 
variations 

5r, (SsxT, Wsx')"- (53) 
With this choice of independent variations, ( ff9| ) yields the conditions: 

- vanishing the individual configurational tractions (|42"|); 

- continuity of the normal Cauchy traction (Eq). 



Additionally, the continuity of the grand canonical potential (HQ) is still equiva- 



lent to the continuity of the normal Cauchy traction (|45|) , by relations ( 42 ) and 
©■ 

The relations ([l2|) , (|4^) for incoherent interfaces with the class of holonom- 
ically constrained variations correspond precisely to the equilibrium results of 
Cermelli & Gurtin (1994) in the spatial description. The equivalent relations 
Iff), M2) correspond to the expressions derived by Larche & Cahn (1978). 



• Coherent interfacial conditions: There are no surface dislocations associ- 
ated with a coherent interface. Furthermore, variations that keep the interface 
coherent should preserve the continuity of x a t S, so that we restrict the class 
of variations to those such that 

IW1 = 0. (54) 
Thus the independent variations at S are 

Sr, 5 sX l - (55) 
With this choice of independent variations, ( f49|) yields the necessary conditions: 

- continuity of the configurational traction 

[c,-]-n s =0; (56) 

- continuity of the normal Cauchy traction (|45|). 
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These two conditions are equivalent to 



- continuity of the Cauchy traction 

le* ® t(]n s = 0; (57) 

- continuity of the normal configurational traction^ 

n s -lE i ®c i K = [c i ]-n 5 =0. (58) 

In this case, the continuity of the grand canonical potential does not generally 
hold. The independent relations ( |5^ ) and ([45]) (or, equivalently, (|57]) and (|58|)) 
correspond to those derived by Cermelli & Gurtin (1994) in the context of 
nonlinear elasticity. These results are also equivalent to the more standard 
expressions based on a fixed reference configuration as originally derived by 
Eshelby (1970) and later by Robin (1974), Larche & Cahn (1978), and GrinfePd 
(1981). 

4.2 Crystal-melt equilibrium 

We now turn to the conditions for crystal-melt phase equilibrium. We identify 
the a-phase as the crystalline solid and denote the melt quantities with the 
superscript F. 

The external boundaries are assumed impermeable and held fixed. Thus the 
total mass of the system is constant: 

M = M a +M F = / ngdv+ p F dv, (59) 
JB" Jb f 

with p F the density per unit volume of the fluid. We assume that the total 
internal energy of the system is given by 

E = E a +E F =( nedv+ [ e F dv, (60) 

JB° JB F 

where e F is the internal energy density per unit volume of fluid, and any surface 
energy has been assumed negligible. Furthermore, we assume that the entropy 
is additive: 

S = S a + S F = [ nsdv+ [ s F dv, (61) 
ie° Jb f 

where s F is the entropy density per unit volume of the fluid. Our constitutive 
assumptions are that: 

e = e a (e i ,g i \ S , Q ), t F — e F (s F ,p F ). (62) 

As before, our condition for equilibrium is 

SQ = 5E -6 B 5S - p 8M = 0. (63) 

8 The normal configurational traction is also called the driving traction at the interface. 
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4.2.1 Unconstrained variations 



• Bulk conditions: We take 

S P F , Ss F , Sg, 8s, Se l (64) 

as the independent variations in the corresponding continuous bulk regions. 
With this choice of independent variations, the analog of ( ^9| ) yields the follow- 
ing independent conditions: 



uniform chemical potential 



n ^ = ^ = ^ (65) 



uniform temperature 



de a de F 

n-^r- = = 0b; (66) 
os os 

dislocation couple balance in the crystalline solid phase 

c 4 +curl k, = 0. (67) 



As before, the couple balance ( p7\ ) implies immediately the configurational force 
balance in the crystalline solid phase which itself is equivalent to the Cauchy 
force balance in the crystalline solid phase. But if we restrict the admissible 
lattice variations to the holonomically constrained class, then configurational 
force balance or Cauchy force balance replaces the couple balance ( |67j ) as the 
mechanical equilibrium condition. 

• Interfacial conditions: We choose 

8r, 5 ei (68) 

as the independent variations at the interface S. They allow the creation and 
transport of defects in the crystalline solid phase. We obtain: 

- continuity of the grand canonical potential 

mo = lu f ; (69) 

- vanishing of the tangential couples of the crystalline solid 

Pk, = 0. (70) 
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Additionally we obtain as a consequence of QTC 
tional traction of the crystalline solid 

c, • n s =0. 

The continuity of the normal Cauchy traction 
rig • (e l ® tj)n s = u F = e F - 



the vanishing of the configura- 

(71) 



0BS 1 



fig 



(72) 



is equivalent to ( |69| ) at S provided (70) and the Eshelby relation ( |3l| ) for Cj. 

Again, if we restrict the admissible lattice variations to the holonomically 
constrained class, then (|7l|) replaces ( |70|) as the equilibrium conditions at S. In 
this case, the relations ([39|), ( |7l|) correspond to those derived by Larche & Cahn 
(1978). 



5 Equilibrium Equations with Interfacial Energy 

The presence of dislocations at the interface as well as the crystalline structure of 
the material may modify the equilibrium conditions at the interface. One way to 
take this into account is to introduce an interfacial energy density that depends 
upon the local state of the interface. We now extend the previous results to 
include such an interfacial energy density. Since the conditions for the bulk 
region do not change, we consider only the appropriate interfacial conditions. 

5.1 The interfacial energy density 

One possible choice of variables to describe the local state of the interface is 

(eT, (eY, n s , (73) 

which reflect the local lattice structure on either side of the interface and the 
orientation of the interface with respect to the current lattice. When the Born 
rule holds, this choice becomes 

(grad X )«, (gradx)' 3 , n s , (74) 

which corresponds to the variables used by Cermelli & Gurtin (1994). Unfortu- 
nately, the dependence on the actual orientation is misleading, since n s does not 
properly account for the crystallographic orientations of the individual lattices 
relative to the interface. A choice of variables that better captures the relevant 
physics at the interface is 



where 



Pie*], n Q , (75) 



K,.( ei ) a E>V |n a .( ei )0Ei 
The choice of variables (J75|) is motivated by the observation that 
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• P^e 1 ]] is the dislocation content of the interface (cf. (fiOj)); 

• n" and n' 3 are the unit normals to the interface in the reference lattices 
of each phase. While it is standard in theories of coherent interfaces to 
include the contribution of the orientation of the interface normal in the 
superficial energy, for incoherent interfaces the energy depends on the 
orientation of the interface normal with respect to both lattices; 

• the fields (P[e 4 ]], n a , n' 3 ) on S may be assigned independently from each 
other. 

Again, the Born rule simplifies ( f75| ) to 

PJgradx], n Q , (77) 



where now 



((gradx) T ) Q n s p _ ((gradx) T ) P n s 



n ' = , ' a - ■ (78) 



|((grad X )-T)«„ 5 |' | ((gradx 

The variables j77| ) correspond to a special case of those used by Leo and Sekerka 
(1989). 

Henceforth, we shall restrict the energy density function to have the special 
form 

^ = p(P|[e < ],n« n"). (79) 

For calculations, however, we will find it more convenient to use the more general 
energy 

i P = t p((e i r,(e i ) ,n s ) (80) 



and then insert ( |79[ ) at the end of the calculation. For example, we denote by 

-- dr " «f-A <-#-. ( 81 ) 



the derivatives of ( j30"l ) with respect to its arguments, which are essentially surface 
stresses. They must however satisfy the following relations obtained by the chain 
rule and ©: 

d<P n s • ( ei ) a ( d(p k 



K a — _P _L £ L_L pa . ( ei ,) a 



k p = P li s ^ % < pP LLL . E fc fei r 

1 0(P|[e*]) |n s .( ej yiV| i ^ J 1 feJ ' 

« - ( po ^ ' E ") M " + ra^ ( p0 ^ ' E *) ^ 

-(°-' e 'l) p ^-(aMll)■ n P|o ■ , • (82) 
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with projections 

P Q :=l-n a ®n Q , P p := 1 - n p ® n 13 . (83) 

An immediate consequence of ( jB2|) is that the surface stresses are tangential 
to S: 

Kf • n s = «f • n s = 0. (84) 
5.2 The variational calculation 

When interfacial structure is taken into account, the grand canonical potential 
( p6| ) includes a contribution from surface energy: 

Q := E + I -6 B S - fiM with I = [ tp da. (85) 

Js 

The variation of the interfacial energy is given by 

5 j <pda = J {< • (fee*)" + «f • (fee*)" + (div s £ - ^) <5r} da, (86) 

where <5s is the variation following the interface (i.e., 8s f = 8f + grad ns / = 
8f + (grad/)n s for any bulk function restricted to S), K is the mean curvature 
of S, grad s and divs are the superficial gradient and divergence on S, and we 
have used the identities (cf. Leo and Sekerka, 1989) 

8 / tpda = I (Ss^p — KipSr) da, #sn 5 = — grad s Sr. (87) 
Js Js 

Remember that the interfacial stresses Kf , and £ are required to satisfy (|82|) . 



With (86) added to the surface integral in (p2j), the interfacial terms in the 



da. 



variation of the grand canonical potential ( |85| ) can be expressed as 

J { (divs £ - Kip - {nujj - [n s x k; • grad ns e 1 }) Sr 

+ < • (S s eT + k? ■ (fee*)" - [n s x k, ■ (fee*)]} 

5.2.1 Case I: unconstrained variations 

When variations are unconstrained in the sense of Section 4.1.1, dislocations are 
allowed to be created and to move both in the bulk and at the interface. We 
choose as the admissible class of independent variations at such an incoherent 
interface 

Sr, (S s eT, (Sse l f. (89) 

This set is more convenient than (|39| ) due to the additional interfacial structure. 
The interfacial equilibrium conditions which replace (40) and ( ^l|) are thus 
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- interfacial normal force balance 

div 5 i-Kip- H + < • (grad^e 1 )" + «f • (grad^e^ - 0; (90) 

- interfacial couple balances 

«f - n s X (k 4 ) Q = and /cf + n s x (k^ = 0. (91) 

The interfacial relations ( |9ll ) state that the couples n s x (kj)"'^ acting on inter- 
facial dislocations and due to the short-range interactions with bulk dislocations 
must balance pointwise the forces n"' 13 in the interface due to interfacial dislo- 
cations. 

The surface divergence of (^) along with the symmetry of the curvature 
tensor imply the 

- interfacial configurational force balances 

di V< s<-(c,) Q -n s =0, and div 5 «f + (c 4 )' 3 • n s = 0, (92) 

which generalizes ([42]), However these are not independent conditions, but 
follow as simple consequences of the interfacial couple balances (^lj). Relations 
( p2[ ) are force balances that state that the "elastic" configurational traction 
acting on any portion TZ of the interface from each phase must be balanced by the 
the interfacial contact force, which tends to reorganize interfacial dislocations 
in 1Z according to their mutual interaction. 

5.2.2 Case II: holonomic variations. 

In this case the variations Sse 1 may be expressed in terms of the inverse defor- 
mation through (|46|) and, by the identity 

P<W = grad 5 (^ X l ) - (grad 5 (n s Sr)) T e\ 

the surface integral in the variation of the grand canonical potential ( p5| ) be- 
comes, granted (§3), 

[div 5 K f - (c t ) a ■ n s ] (S SX T + [div 5 «f + (cif ■ n s ] (S sX Y 
- (div S £-K<p + (e l )« • n s div s < + (e l f • n s div 5 «f 

+ k? • (grad ns e J )« + «f • (grad^e*)" - n s • [e* ® t 4 fln s ) Sr] 0$) 
where Cj is now given by (j50j). 

• Incoherent interfacial conditions: We take 

Sr, (S s xT, MY (94) 

as the admissible class of independent variations. This yields as generalizations 
of ( f42| ) and (^3) the interfacial equilibrium conditions: 
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- interfacial normal force balance 

div 5 £ - Ktp + (e l ) a ■ n s div 5 < + (e 1 )' 3 • n 5 div 5 »f 

+ < • (grad ns e 4 ) Q + (k^ ■ (grad ns e^ - n s • [e* ® t,]n fl (95) 

- interfacial configurational force balances ( |92| ) . 

Now however (|9^ ) are independent conditions. Furthermore, ( |95|) is equivalent, 
granted ©, to ©. The equilibr ium interfacial conditions (|90|) and (|92|) arc 
appropriate to a constrained situation in which dislocations cannot be created 
or move in the bulk and, for this case, arc equivalent to the relations derived by 
Cermelli & Gurtin (1994). 

The non-holonomic equilibrium conditions ([Jl]) are stricter than their holo- 
nomic counterparts (|9^); interfaces which are in equilibrium governed by (|92|), 
i.e., when dislocations are absent in bulk, need not be in equilibrium if disloca- 
tions were allowed, since these dislocations could exert forces on the interface. 
As an example, consider a grain boundary between undeformcd grains in a 
polycrystalline material, and assume that the grains are rotated one relative to 
another: 

(e,)" = Ei, (e,)' 3 = QE„ (96) 

with Q a rigid body rotation. Assume that dislocations are not allowed in bulk 
(holonomy), fix g — go and s = so, and take the grand canonical potential 
such that, by objectivity, u)(Ei, go, sq) = tj(QEi, go, sq) = is a minimum, so 
that { Ci ) a = ( Ci f = (cf. @). Then, according to the holonomic equilib- 
rium condition (|9^), any interface such that /t" ,/3 is a constant (for instance 
planar interfaces with a uniform dislocation density) is in equilibrium. On the 
other hand, if dislocations are free to move to and from the interface from the 
bulk, equilibrium is established by the non-holonomic relations (|9l|), and such 
an interface (planar with a uniform dislocation density) may well not be at equi- 
librium, as for instance when k" is non-zero but kj is zero. 



5.3 The equilibrium conditions projected on a slip system 

The transport of defects is commonly interpreted in terms of a crystallographic 
slip system, where slip typically takes place along lattice vectors. Some insight 
into the structure of the bulk and interfacial equilibrium conditions may be 
gained by introducing such a crystallographic slip system and projecting the 
equations onto this slip system. The resulting relations involve forces acting 
within the slip plane, which may be interpreted as forces on the dislocations 
moving by glide on such a plane. 
More precisely, let 

tr<8>fi, (97) 
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be a slip system, with cr a slip vector and \x the slip-plane normal, with the 
restriction that fi ■ cr = 0. The slip vectors are taken to be lattice vectors, which 
correspond to the Burgers vectors of the dislocations whose motion generates 
the slip, while the slip-plane normals are reciprocal lattice vectors. 

The slip system is measured in the reference lattice. Its counterpart in the 
current lattice is represented by the set of vectors 

cr'm, with cr i :=<r-E i , (98) 

and 

m = fj,je > , with fij := /x • R,-. (99) 

Note that a 1 and fj,i are constants. 

Projecting the bulk dislocation couple balance ( |36| ) onto the slip system (Pq), 
we obtain 

a' Ci ■ m + ctVj (ki ■ g j ) + a' div (kj x m) = 0, (100) 



with the configurational stress given by (j3l|). Each term in ( 100 ) has a suggestive 
interpretation, namely: 

• the term a 1 C; • m is the configurational shear across the slip plane (recall 
in fact that m is the current slip-plane normal, while cr is contained in 
the referential slip plane) , and may be viewed as forcing or inhibiting slip 
along the corresponding system. It may be identified with the so-called 
long-range stress on dislocations in a dislocated crystal. 

• the term a % x m is the couple acting in the slip plane on dislocations with 
Burgers vectors parallel to cr. It may be identified with the short-range 
contact couple due to the presence of other dislocations. 



Projection of the interfacial couple balances (91) onto a phase-a slip system 
yields 

(ctT m • < - (ctT ((ki) a x m) • n s = 0. (101) 

The couple m x (k^)" is tangential to the slip plane, while Kf is tangential to 
the interface. Consider now the intersection line (say L) of the slip plane with 
the tangent plane to the interface. A straightforward calculation shows that 

• the term (a l ) a m • nf represents the projection onto the direction m — 
(m ■ n s )n s , the perpendicular to L in the tangent plane to the interface, 
of the interfacial force on a dislocation with Burgers vector cr. 

• the term (a l ) a (m x (k^) a ) -n 5 represents the projection onto the direction 
n s — (m • n s )m, the perpendicular to L in the slip plane, of the couple 
due to bulk dislocations on a defect with Burgers vector cr. 
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A dislocation lying along the intersection of the slip plane and the interface is 
parallel to L\ the only relevant component of any force acting on this dislocation 
is the component orthogonal to L. The couple balance requires that these 
components of the opposing couples on the dislocation be equal. Therefore ( 101 ) 
represents a balance of couples which tend to redistribute the dislocation density 
at the interface and couples which tend to move the dislocations away from the 
interface along their slip planes, due to the presence of bulk dislocations. Two 
such independent balances are needed, one for each phase, since dislocations 
may be transferred to and from the interface from each crystal phase, and their 
slip systems are independent. 



6 Incremental Deformations (without Interfacial 
Energy) 

In many models of defective materials, it is common to introduce an elastic- 
plastic decomposition of the macroscopic deformation. In our formulation there 
is no need to introduce the notion of deformation from a fixed reference con- 
figuration. In order to compare the two approaches, however, we now intro- 
duce such an elastic-plastic decomposition. We follow the procedure of Davini 
(1989) by superposing a conventional elastic deformation onto a dislocated crys- 
talline solid, where the elastic deformation does not change the detectivity of the 
crystal. We then derive the necessary equilibrium conditions for the resulting 
incremental theory. 

We identify B, which corresponds to the region occupied by the crystalline 
solid in a dislocated state, to the intermediate configuration in an elastic-plastic 
decomposition. In contrast to the standard formulations of the elastic-plastic 
decomposition, we do not assume that B is stress-free. Now consider a defor- 
mation f : B — > M 3 from this dislocated state that is continuous but possibly 
non-smooth across the surface <S. Although the incremental deformation is con- 
tinuous everywhere, the interface can be incoherent due to the presence of sur- 
face dislocations in the dislocated intermediate configuration. We assume that 
the lattice vectors transform according to the Born rule under the deformation 
f: 

e 4 i ► e ; := Fe l5 e i i-> e l := F~ T e\ (102) 

where the Sj are the current lattice vectors in the deformed region f(£>) and 
F = grad f . The dislocation density associated to the field e l is given by 

g i = det(F- 1 )Fg i , (103) 

which yields 

~r U>- (104) 

Thus the dislocation density per unit cell is unaffected by the incremental elastic 
deformation. 
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The lattice variations of this resulting state can be expressed as (Se^Sg), 
which, as before, are assumed to vanish on the external boundary dB. Consider 
now the additive composition of the two configurations 



-i A 



XSe u 



Qx := g + XSg 



(105) 



where A is a small parameter. This composition induces a corresponding varia- 
tion in any functional $ depending on the composed configuration: 



d 

5 * = aA*[ eiA ' 0A - 



A=0 



Note however that 



<5e, = (SF) ei - (e; • *e*)Ete* 



(106) 



(107) 



Because of (102) and (104), we may assume that the energy densities are 
functions of the form 



e = e a (ei(ei, F), g iJ , g, s) for phase a, 
e = e^e^e^F),^- 7 , g, s) for phase (3, 



(108) 



where the unbarred quantities refer to the dislocated configuration B. We in- 
troduce the incremental Piola stress 



Tr := n—r (g) e; 
oei 



de_ 



so that the Eshelby relation (31) becomes 
d(nui) 



de l 



no; e, - T R Fe 4 - (g J x e,) x k, 



(109) 



(110) 



due to the superposed deformation. 

Thus the expression ( p8|) for the variation of the energy functional becomes 



Sfl 



[n{^- - fi)Sg + n(^- - 6 B ) 5s 



+ (ci + curl ki) ■ Se l - (div T R ) • Sf] dv 

- [ (ImJ - n s • F T T R nJ Sr + n s ■ {Se* x k<] 
Js 

+ [T R nJ-fef) da = 0, 
where we have used the identity 

-lT a n s • Sf} = -[T R n s ] • 6 s f + n s ■ [F T T R nJ Sr, 
and introduced 

S s f := Sf + (grad ns f) Sr = Si + Fn s Sr, 



(111) 



(112) 



(113) 
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which is the variation of the deformation f following the variation of the inter- 
face. Notice that the continuity off implies the continuity of 5$f. The difference 



between (111) and ( p2| ) is due to the superposed elastic deformation f. 
We take as the admissible class of independent variations in bulk: 

5q, 5s, 5f, Se\ (114) 

The introduction of the incremental deformation f provides an additional ad- 



missible variation to (|33|). This choice in (111) leads to the following necessary 
conditions :[] 



uniform chemical potential 



de a _ del 3 
dg dg 



uniform temperature 



3e a de 13 



ds ds 
Piola force balance for each phase 

divT R = 0; 

dislocation couple balance for each phase 

C; + curl = 0. 



(115) 



(116) 



(117) 



(118) 



As in the previous unconstrained cases, the couple balance (118) implies imme- 
diately the configurational force balance for each phase 



div c, = 0, 



(119) 



but, importantly, it is not equivalent to the force balance ( 117 ). The inde 



pendent conditions in the bulk are (fTTq)-(llS). In particular, there are two 



distinct mechanical equilibrium conditions: the couple balance (118) enforces 



mechanical equilibrium of the dislocated intermediate region 6; the force bal- 



ance (117) enforces mechanical equilibrium of the elastically deformed region 
f($). Interestingly, both the intermediate and the current configurations must 
be in mechanical equilibrium. p| 

9 The derivatives are, as before, with respect to the position x, which is now in the reference 
configuratio n fo r the i ncrem ental elastic deformation. 



10 Both (117) and (118) correspond to equations derived variationally by Le & Stumpf 
(1996) in a continuum theory of dislocations with a postulated elastic-plastic decomposition, 
though they were not interpreted as the mechanical equilibrium conditions in the actual and 
intermediate configurations. 
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If we restrict the admissible lattice variations to the holonomically con- 
strained class (i.e., the Born rule) then (119) replaces (118) as the mechanical 
equilibrium condition in the intermediate configuration. Then both the config- 
urational force balance (119) and the standard force balance (117) are needed 
to enforce mechanical equilibrium. 

At the interface we take as the admissible independent variations 



6r, S s f, (Se l ) a , (5e^ . 
They yield the following interface conditions: 

- continuity of incremental normal Eshelby traction 

n s ■ Inul - F T T R ]n s = 0; 

- continuity of incremental Piola traction 

[T R ]n 5 = 0; 

- vanishing of the individual tangential couples 

P(k,) Q = and P(k l ) /3 = 



(120) 



(121) 



(122) 



(123) 



Both (121 ) and (122) are standard equilibrium conditions for coherent interfaces, 
whereas (123) is the condition for mechanical equilibrium for an incoherent in- 
terface. Both sets of conditions arise because the interface is coherent under the 
superposed elastic deformation, but incoherent in the intermediate dislocated 
configuration. As in the bulk, we see that the interface must be in equilibrium 
in both configurations when an elastic-plastic decomposition is introduced. 

If we restrict the lattice variations in the intermediate configuration to the 
holonomically constrained class, then the vanishing of the individual configura- 
tional tractions in the intermediate configuration replaces the equilibrium con- 
dition (123). Note that this configurational stress differs from the incremental 
Eshelby stress, which arises from the elastic part of the deformation. 



7 Discussion 

One way of introducing the notion of a discrete crystalline lattice into a macro- 
scopic continuum model is through additional fields that represent, for example, 
the lattice basis vectors of the crystalline structure at each point. It is com- 
mon however to eliminate these additional degrees of freedom by means of the 
Born rule, which states that the lattice vectors at each point deform with the 
gradient of the macroscopic deformation. Thus a crystalline body is reduced 
to an anisotropic elastic continuum, which remembers its lattice structure only 
through the symmetry of the constitutive relations. This approach fails when 
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defects such as dislocations have to be accounted for in the model, since the 
Born-rule no longer holds. In this case the macroscopic deformation gradient 
fails to describe the true local ordering of the lattice. 

In this work, we use the kinematic ideas of Davini and Parry on defective 
crystals to construct a model of crystalline solid equilibrium. The basic assump- 
tion is that the defect density is sufficiently small, so that an average ordering 
is still recognizable locally. At each spatial point of the solid a lattice basis may 
therefore be assigned, but these local descriptions need not match globally to 
yield a smooth deformation of a perfect crystal. Dislocations are then viewed as 
the obstruction to patching together such local lattice bases to a global lattice 
basis. 

Using a variational framework, we obtain equilibrium equations for the bulk 
and for incoherent interfaces in the presence of vacancies and dislocations with- 
out assuming the validity of the Born rule. Although the mechanical equilib- 
rium equations involve dislocation couples, indeed they also imply the more 
familiar equilibrium conditions ensuing from the Cauchy and configurational 
force balances. We discuss how these equations simplify when the Born rule 
is imposed as a constraint on the possible lattice variations so that the defect 
density vanishes in bulk, and also when coherent boundaries are considered. In 
these constrained cases, our results correspond to those previously derived by 
Larche & Cahn (1978), Leo & Sekerka (1989), and Ccrmelli & Gurtin (1994). 
Importantly, our unconstrained results provide a generalization to their work, 
which only dealt with conventional elastic bodies, or, in our language, holonomic 
microstructures. 

In order to compare to common models involving an elastic-plastic decom- 
position of the strain, we also discuss an extended model, in which we superpose 
conventional elastic deformations onto equilibrium states of a dislocated crystal. 
The resulting system for the bulk is composed of a force balance that enforces 
mechanical equilibrium of the current configuration, together with a couple bal- 
ance enforcing the mechanical equilibrium of the dislocated intermediate con- 
figuration. Such bulk conditions have also been derived by Le & Stumpf (1996) 
in a slightly different context. The interfacial conditions are composed of the 
standard coherency relations requiring the continuity of the Piola traction and 
the continuity of the usual Eshelby-stress traction, but complemented by an 
incoherent relation requiring the vanishing of the individual tangential couples, 
which enforces the equilibrium of the interface with respect to the lattices of 
each phase in the intermediate dislocated configuration. 
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